Unconventional photon blockade in coupled optomechanical systems 
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We show that in laser-driven coupled optomechanical systems, photon antibunching can occur un- 
der weak optomechanical coupling, contrarily to common expectation. This unconventional photon 
blockade originates from destructive quantum interference between different excitation pathways 
bringing from the ground to two-photon states. Using a quantum open-system approach, we study 
the antibunching as a function of driving field amplitude, temperature, and pure dephasing rate, 
and derive optimal values of the system parameters for its occurrence. These values are remarkably 
close to those characterizing optomechanical systems in some current experimental studies. 

PACS numbers: 42.50.Wk, 42.50.Ar, 03.65.Yz, 42.65.-k 



Optomechanics is the study of optical cavities in which 
the resonant mode is coupled to a mechanical oscilla- 
tor by means of radiation pressure. The field has ex- 
perienced a tremendous progress in recent years [U-OI], 
marked by several achievements, among which the suc- 
cessful optical cooling close to the ground state of the 
mechanical motion [J-lff and the coherent coupling be- 
tween the two degrees of freedom . In addition to the 
promise of operating optomechanical systems (OMSs) as 
extremely sensitive detectors of mass, force or position, 
perhaps the most interesting perspective is the realiza- 
tion of strong single-photon nonlinearities. These would 
bring to several possible schemes for producing quantum 
correlated states of light and macroscopic mechanical de- 
grees of freedom, with possible applications in quantum 
information processing 0, |TlJ . 

An OMS is predicted to display single-photon non- 
linearities |1-12|, provided the optomechanical coupling 
is strong enough to overcome system losses - in anal- 
ogy to what hap pens in cavity quantum electrodynam- 
ics (CQED) fl3|. In particular, in both single El flT 



and coupled [ljj, LL2| OMSs, the emission of photons with 
sub-poissonian statistics is expected. This kind of photon 
blockade is what occurs for a resonantly driven two-level 
system after emission of one photon, the system 

needs to be excited again and a second photon emission 
is most likely to be delayed by one lifetime. Any system 
with a strongly anharmonic set of energy levels, when 
driven resonantly to one of its possible optical transi- 
tions, will essentially behave in the same fashion. 

A mechanism of completely different nature however - 
that has been recently named (ljj unconventional photon 
blockade (UPB) as opposed to the conventional mech- 
anism (CPB) described above - can give rise to sub- 
poissonian photon emission in CQED even in presence 
of much weaker nonlinearities. In this mechanism, the 
presence of two photons in the system is prevented by 
destructive quantum interference - enforced by the small 
nonlinearity - between different excitation pathways that 
can lead from the ground to the two-photon state. UPB 
has been initially predicted by Carmichael [HI for the 



Jaynes-Cummings model and successfully demonstrated 
in atomic CQED experiments fl7l Il8j |. More recently, 
UPB was predicted 3, 2{| in coupled CQED systems in 
presence of very small nonlinearities - both of the Kerr 
and Jaynes-Cummings types - such that the phenomenon 
might even be induced by the background material non- 
linearity of modern semiconductor nanocavities plj . 

In this letter, we investigate the occurrence of UPB 
in coupled OMSs, and characterize it as a function of 
losses, driving field strength and dephasing mechanisms. 
We show that UPB is expected to occur well in the weak 
optomechanical coupling limit, contrarily to CPB. By re- 
viewing the state-of-the-art parameters of current OMSs, 
we show that some of them are close to the UPB re- 
quirements and set up a perspective for the experimental 
observation of this phenomenon. 

Consider two linearly coupled optical modes, described 
by bosonic creation operators a\ and a\. Mode two is 
characterized by optomechanical coupling to a phonon 
mode described by a bosonic creation operator b\. Mode 
one is driven by a monochromatic laser field at frequency 
Lu p . In the frame rotating with the driving frequency, the 
Hamiltonian is (assuming h = 1) 



H 



Ab 2 



Aiajai + A2S2S2 - 
J(a\a 2 + a\ai) + ga\a 2 (b\ + b 2 ) + e(a\ + &i)(l) 



where Ai = lo\ — lu p and A2 = co 2 — lo p are the detunings 
of the two mode frequencies, J is the rate of the linear 
coupling, uj m is the mechanical frequency, g is the single- 
photon optomechanical coupling rate, and e denotes the 
(real) amplitude of the driving field. We further assume 
dissipation rates K12 and 7 for the optical and mechanical 
modes respectively. 

The dynamics of this driven-dissipative system is gov- 
erned by the Von-Neumann equation for the density ma- 
trix, p(t): 



dp 
If 



3=1,2 

(N th + l)V[b 2 ]p + N th V[bl]p 



(2) 




where V[6]p = 20pO t - O^Op - pO^O are Lindblad 
terms describing the markovian interaction with the en- 
vironment and Nth = (e Um / KBT — is the average 
occupation number of the mechanical mode for a given 
temperature T. 

In the absence of driving terms, the photon and phonon 
degrees of freedom in an OMS can be decoupled by 
means of a polaron transformation (9[. As a result, 
the optomechanical coupling term in ([I]) is replaced by 
H' = —A g (a,2a2 +02026202), with A g = g 2 /uj m , repre- 
senting a polaron energy shift and an effective Kerr-type 
nonlincarity on mode two. This sets a close analogy be- 
tween the present system and the weakly nonlinear cou- 
pled optical modes studied in Refs [lilU, that were 
shown to display UPB under appropriate conditions on 
the system parameters. We thus set out to demonstrate 
that Hamiltonian (JTJ can also give rise to UPB and in- 
vestigate its features. 

We assume Ki = K% = k, and use this quantity to 
rcscale all other energy parameters. The UPB arises 
in presence of a weak nonlincarity in the second opti- 
cal mode, which amounts here to choosing A 9 /k <C 1. 
It further requires u) m / 'k S> 1 - namely the well-resolved 
sideband limit - as explained below. These two con- 
ditions naturally allow for a weak optomechanical cou- 
pling g/n < 1. We consider three fixed values of the 
mechanical frequency: uj m = 11k, 24k, and 50k, charac- 
terizing the state of the art in different kinds of OMSs 
0, El EH- UPB depends critically on the value of A2 



(while being much less sensitive to Ai [20t l21|). Hence, 



throughout this work, we set the value of the polaron- 



shiftcd detuning A2 
tion derived in Rcf. 



- A g to the optimal UPB condi- 
[13 in the e -> limit, A 



Opt 



-(1/2)vV9J 2 + 8k 2 J 2 



3J 2 . We additionally set 



Ai = A opt and 7 = 0.01k, and study the result as a 
function of the remaining parameters. 

For the numerical analysis, we assume steady state, 
p = 0, and solve the linear system arising from ^ un- 
der the condition Tr(p) = 1, by truncating the Hilbert 
space of the system to the basis 712, n m )} with 

n\ + tt-2 < N p h and n m < N m . The numerical solu- 
tion is obtained by means of a highly-optimized numerical 
code (see supplemental material). We have numerically 
checked that our results are scarcely affected by the as- 
sumption of an additional mechanical mode b\ coupled 
to the first optical mode. This is expected, as the pho- 
ton occupation is much larger in mode two 19|, 2^. For 
all results, convergence vs N p h and N m was accurately 
verified. 

In the present scheme, we expect the occurrence of pho- 
ton antibunching in the first mode. We therefore com- 
pute the quantity g( 2 '(0) = (a\a\aiai) / (a\di) 2 , where 
the quantum expectation value is evaluated as (O) = 
Tr(pO). Fig. ma) is a color plot of log 10 ( 5 ( 2 ' (0)) as a 
function of g and J, for uj m = 50k and a weak probe 




FIG. 1: (a) log 10 (s (2) (0)) computed for OJm = 50k and e = 
0.1k, as a function of g and J. White line: optimal UPB 
condition given by Eq. (5) of Ref. (b) Three horizontal 
cuts of the color-plot in (a) for the values of J indicated in the 
legend. Dashed lines: same quantity computed by restricting 
to the zero-phonon subspace. 



amplitude e = 0.1k. For most values of the coupling rate 
J, antibunching occurs in a well defined range of values 
of the optomechanical coupling (blue in the plot), with 
g^ 2 \0) taking values as small as 10~ 3 . The white curve 
in the plot represents the optimal UPB condition on A ff 
derived by Bamba et al. for the case with Kerr nonlin- 
earity [19( (see Eq. (5)). The agreement with our result 
confirms that the antibunching is indeed produced by the 
same UPB mechanism. Fig. [IJb) shows g^ 2 '(0) as a func- 
tion of g for three selected values of J. We observe that 
a significant antibunching can be achieved for values of 
the optomechanical coupling as small as g = 0.1k. 

The present result differs however substantially from 
the one with Kerr nonlincarity in that the minimal value 
of 5 ( 2 )(0) for a given J is subject to a lower bound. In 
the Kerr case, for the optimal choice of the parameters, 
<r 2 )(0) decreases with the intensity of the driving field 



and can become arbitrarily small 3, 1^, 2(|. To under- 
stand this behaviour, we study the limit of vanishing e by 
restricting the space of relevant Fock states \ni,n 2 ,n m ) 
to Nph = 2 and N m = 1 (implying a total of 12 states). 
The corresponding bare energy levels are sketched in Fig. 
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FIG. 2: Sketch of the energy levels corresponding to Fock 
states up to n p h = 2 and n m = 1. Black (red) horizontal 
lines denote zero (one) phonon states. Dashed arrows indicate 
possible excitation pathways to states |200) (black and green) 
and |201> (red). 




m 



FIG. 3: Thick lines: minimal value of <? (0) obtained by 
varying g, plotted as a function of J. Thin lines: correspond- 
ing value of g for which the minimum occurs. 



[2]together with the relevant coupling terms. In order to 
achieve UPB, the amplitude of both the |200) and |201) 
states must be suppressed. In the case with Kerr non- 
linearity, the |200)-amplitude is suppressed via destruc- 
tive interference between the linear excitation pathway 
(sketched as black dashed arrows) and all those involving 
the nonlinearity arising from an even number of phonon 
creation/destruction processes (e.g. the green pathway 
in the sketch). This condition is achieved for the opti- 
mal conditions derived in Ref. [l9|. In the present case 
however, there are also excitation pathways leading to 
the 1 201) state, involving an odd number of phonon cre- 
ation/destruction processes. The same optimal condi- 
tions cannot suppress both amplitudes simultaneously, 
which ultimately explains the lower bound on g^ 2 \0) 
observed in the numerical results. The existence of 
this lower bound can be proved analytically by studying 
the Schrodingcr equation on the 12-state basis described 
above (see supplemental material). To further support 
this picture, using the full model ([2]), we evaluate the sec- 
ond order photon correlation under the condition of zero 
phonon occupation. This can be defined using the oper- 
ator Po = I\ ® I2 ® |0 m )(0 m |, that projects on the zero- 
phonon subspace: <?q 2 ^(0) = (a\a\aiaiP ) / (a\aiPo) 2 . 
For each of the three sets of parameters considered in 
Fig. [TJb), this quantity is plotted as dashed red lines. 

The lower bound on g^ 2 \0) is clearly suppressed and the 
(2) 

minimal value of <?q (0) is now set only by the finite value 
of the driving field e used in the numerical calculation. 

In order to quantify the effectiveness of the UPB mech- 
anism, we have minimized the value of g^ 2 \0) by vary- 
ing the optomechanical coupling g. Fig. [3] displays 
min(g( 2 )(0)) (thick lines), and the corresponding value 
of g (thin lines), for varying J and for the three values 




FIG. 4: g (2) (0) as a function of the average photon occupa- 
tion of the first mode m. Circles: parameters g and J set 
to the absolute minimum of g (2) (0) in Fig. [3] Tr iangles: pa- 
rameters g and J corresponding to min(g^ 2 ' (0)) = 0.1 in Fig. 
[3] (points corresponding to different values of u) m sit on the 
same curve). The rightmost points required up to N p h = 22 
for convergence. 



of u) m considered. The absolute minimum of g^ 2 \0) is 
reached respectively for J = 1.9ft (u> m = 11k), J = 2.6k 
{uj m = 24k), and J = 3.6k (u m = 50k), always signifi- 
cantly smaller than cj m , with the corresponding optimal 
value of g ~ 1.1k for the three cases. A small feature 
at 2J = uj m — A2 arises because of resonant enhance- 



ment of the nonlinearity in coupled OMSs [ll|, LL2J , but is 
scarcely significant in the present weak optomechanical 
coupling case. We see that, for the three mechanical fre- 
quencies, g( 2 '(0) ~ 0.1 can be achieved for g ~ 0.1k, cor- 
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responding to A g < 10~ 3 k, and a somewhat larger value 
of the optical coupling J. The increase in mm(g^ (0)) 
with decreasing is a consequence of the interferen- 

tial nature of UPB. In fact, as suggested by the scheme 
proposed in Fig. [2J the well resolved sideband limit corre- 
sponds to having a small relative uncertainty in the phase 
of an excitation pathway involving phonon creation and 
destruction processes, which is in turn necessary for the 
occurrence of destructive interference. 

To assess the expected rate of emission of antibunched 
photons, we study how g^ 2 \0) varies as the driving field 
amplitude e is increased. Fig. [4] displays g^ 2 \0) as a 
function of the average photon number in the first mode 
m = (a\ai). Circles are computed for values of g cor- 
responding to the three absolute minima of g^'(0) in 
Fig. O while squares correspond to values of g and J se- 
lected to have g {2) (0) = 0.1 in Fig. H This latter plot is 
identical for the three values of u) m , showing that </ 2 '(0) 
ultimately depends on the photon occupation only. More 
generally, we observe that g^ 2 \0) preserves its minimal 
value up to n\ ~ 10~ 2 , contrarily to the Kerr case [io| . 
where g^ 2 '(0) decreases as a power law of m. This is due 
to the lower bound set by the phonon-assisted processes 
discussed above. For n\ > 10~ 2 the antibunching is pro- 
gressively suppressed, and the coherent result g^ 2 \0) = 1 
is recovered as the classical limit n\ 1 is approached. 




FIG. 5: 3 (2) (0) vs J for different values of the pure dephasing 
rate V, computed for (a) u m — 11k, (b) u m = 24k, and (c) 
u) m = 50k. For the three cases, g was set to its optimal value 
corresponding to absolute minima of 3 (2) (0) in Fig. El (d) 
(/ 2 '(0) vs temperature, as computed for oj m = 24k and the 
optimal values of g and J. Black line: thermally activated 
behaviour p' 2 '(0) oc exp(A'_BT/2oj m ). The rightmost points 
in (d) required up to N m = 36 for convergence. 

Given the interferential origin of UPB, it is important 
to study how it is affected by pure dephasing induced 
by interactions with the environment. Pure dephasing 



typically arises from scattering processes of one photon 
with emission or absorption of one bath quantum [lij ]. 
For mode one, they are described by the Hamiltonian 
Hp = J2k 9k<i\ai(bk + b\), where Sj, are Bose opera- 
tors describing modes of the thermal bath with given 
frequency uJk- A similar Hamiltonian holds for mode two. 
In the present case, terms of this kind can be interpreted 
as the optomechanical interaction of each photon mode 
with mechanical modes other than the one of interest. 
To lowest order, and within the Markov approximation, 
they are described as additional Lindblad terms to Eq. 
© 



dp 
~dt 



P . 



(3) 



where the magnitude of the rate T - assumed equal for 
the two modes - is set by the polaron energy gl/uik- Fig. 
[5] (a)-(c) show the variation of <r 2 ^(0) as a function of 
J for varying T, for w m = 11k, 24k, and 50k respec- 
tively. As expected for a phenomenon based on quan- 
tum interference, the UPB is suppressed by pure dephas- 
ing. The suppression is complete at a rate T ~ 0.01k, 
namely as soon as T > A g . This is easily understood, 
as A g quantifies the phase shift required for exact de- 
structive interference between excitation pathways to the 
two-photon state. In typical OMSs, the mechanical mode 
under study is tailored to maximize the optomechanical 
coupling g, among all other modes in the same frequency 
range. We thus expect that the condition T <C A g is 
fulfilled in most OMSs currently investigated. 

Finally, Fig. [SJd) shows the variation of g^ 2 \0) with 
temperature, computed for the optimal UPB parameters 
in the case uj m = 24k. The effect of temperature be- 
comes relevant as soon as KqT > oj m , where we see 
the onset of a thermally activated behaviour g^(0) cx 
exp(K B T/2uj m ). Accordingly, we expect g (2) (0) > 0.1 
for KbT > 150k, which is less than one decade below 
the value KbT ~ 1000k characterizing the mechanical 



mode in Ref. 22 



The conditions required to observe UPB are radically 
different from those needed for the convetional mecha- 
nism. Firstly, the monochromatic driving field is not 
resonant with any of the eigenstates of the coupled sys- 
tem (light is absorbed by virtue of the finite linewidths). 
But most importantly, UPB requires a much smaller op- 
tomechanical nonlincarity A s 



io- 



to 10- 



for A 9 /k > 1. 



9ZJu m , ranging from 
In single OMSs Q, CPB is observed 
» s/ n, ^ _l. In coupled OMSs on the other hand, the 
resonance condition u) m ^ 2 J has been shown to enhance 
nonlinear effects also for A 9 /k < 1 [ll|, [l2|. These reso- 
nant schemes however still rely on a strongly anharmonic 
spectrum and are thus bound to the additional condition 
of strong optomechanical coupling g 3> k. Here, we have 
shown that UPB is achieved for the much less stringent 
condition g/n ~ 0.1 -f- 1. These considerations suggest 
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that UPB should be much more easily achievable in state- 
of-the art OMSs than CPB. Candidate systems must be 
in the well resolved sideband limit to m 3> re, thus exclud- 
ing Bose-Einstein condensates of trapped ultracold atoms 
23l l24j for which u m is in the kHz range. To this pur- 
pose, the best strategy is probably to minimize the cavity 
loss rate re rather than to maximize ui m at the expense 
of a lower optical quality factor [iH . Ultimately, the sys- 
tems closest to the UPB requirements are silica toroidal 
microresonators (w m /re = 11, q/re = 5 x 10~ 4 ), op- 

24, 



2J, |26| (w m /re 



tomcchanical crystal nanobeams 
g/re = 5xl0~ 3 ), and micromechanical oscillators coupled 
to microwave cavities [j| (ui m / k = 50, g/n = 10~ 3 ). The 
first two systems easily allow for a wide range of values of 
J/re 2fi|, 27 1 . Optomechanical crystal nanobeams [22| in 
particular would require an improvement of the ratio g/n 
by a factor of 20, and cooling down to T ~ 1 K. Another 
promising system is the Fabry-Perot GaAs/AlAs micro- 
cavity with distributed Bragg reflectors [Hf (w m /re = 1.7, 
g/n = 4 x 10~ 3 , re ~ 27r x 12 GHz), where UPB require- 
ment could be easily reached for the current state-of-the- 
art loss rate re achievable in these systems . 

Quantum interference in coherently driven optical sys- 
tems is a promising mechanism in view of the generation 
of quantum correlated states in presence of weak non- 
linearities 15. 30, [3TI . Applied to OMSs, this paradigm 



might be the shortest route towards the controlled gen- 
eration of quantum correlations between optical and me- 
chanical degrees of freedom. 
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Appendix A: The numerical method 



The numerical calculations in the cases of large driving field amplitude e, or large temperature, required values of 
Nph up to 22 and of N m up to 36 respectively. Solving the stationary von Neumann equation for the density matrix 
p on such a large basis set (up to 800 states), as was done in the present work, requires an optimized numerical 
approach, that we describe below. 

The starting problem is a homogeneous linear system of equations for the elements of the density matrix, 



Cp = 0, 



(A.l) 



subject to the condition Tr(p) = 1. Numerically, we need to "vectorize" the matrix p, namely to express it as a 
column vector vec(/5) made by stacking the columns of p on top of one another. If p is a N x N matrix (we omit the 
hat notation for all matrices below) , then 



vec(p) = [pi 



■ pN,l,Pl,2, ■ ■ ■ PN,2i 



, Pi, JV, ■ ■ • PN,N\ 



(A.2) 



To express the Lindblad superoperator in matrix form, to be used in this representation, we rely strongly on the 
property 



vec (ApB) = (B T <g> A) vec(p) 



(A.3) 



where A, p, and B are Nx N matrices and <£> denotes the Kronecker product. The matrix B T ® A is a N 2 x N 2 matrix. 
Indeed, all terms of the von Neumann equation are expressed as products of operators with the density matrix. In 
particular, the above identity gives 



vec{[H,p}) = (I ® H - H T <g> I) vcc(p) , 
vec(2oj/oaj — a^ajp — pOjdj) = ^2(a]) T <g> aj — I <X> {a\a>j) — (8 1^ vec(p) 



(A.4) 
(A.5) 



where I is the N x N identity matrix. From these identities, the whole Lindblad superoperator C can be cast in 
matrix form. 

To build the matrix corresponding to C acting on a truncated Fock-state space, the following procedure is applied. 
We start from the destruction operators on the truncated Fock-state space of each single mode (the creation operators 
are then simply obtained as the adjoint operators). We have 



/0 1 ••• 

0V2- 



a N ph 



\ 







yiv 





l>N m 



( 1 

y/2 







/AT 




\ 



/ 



(A.6) 



From these matrices we can build the corresponding matrices acting on the truncated Fock-state space of the multi- 
mode system. For example, for the case considered here of two photon and one phonon mode, the matrices of the 
corresponding destruction operators arc built as 



a l = a N ph ® J N ph ® T N m , 

a-2 = lN ph ® a Nph ® lN m , 
b 2 = lN ph ® lN ph ® b Nm , 



(A.7) 
(A.8) 
(A.9) 
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where as before, In is the N x N identity matrix. These matrices are now operating on a (N p h + l) 2 (N m + 1)- 
dimensional vector space. We can further reduce the computational task by restricting to the space generated by 
states |ni, ri2, n m ) such that n\ + n 2 < iVp^ and n m < N m . This step is justified by noticing that n\ + n 2 determines 
the total energy of the photon field, so that the truncation corresponds to retaining all states within a maximum 
total energy, responding to an intuitive perturbation criterion. As an example, for N p h = 22 and N m = 2 we have a 
Fock-state basis of (N p h + l) 2 (N m + 1) = 1587 states, while only 828 states are left after the additional truncation 
procedure. Numerically, this truncation is obtained by building the vector n p h obtained from the diagonal of the 
matrix a\a\ + a3,a 2 , searching all elements of this vector for which n p h(j) > N p ^ and removing from the matrices 
(|A.7[) - (|A.9j) the corresponding j-th rows and columns. 

From the matrices (|A.7[) - (|A.9[) . the Hamiltonian can be computed by simple matrix multiplications, following which, 
the matrix of the Lindblad superoperator C can be computed using the prescriptions (|A.4[) - (|A.5p . Finally, the linear 
set of equations (|A.1| , with the additional condition on the trace of the density matrix, can be solved by replacing 
one of the equations in the set by the condition Tr(p) = 1 expressed in vectorized form. This is possible as we are 
assuming the existence of a solution to the stationary problem, and thus det(£) = 0. The resulting linear system is 
now inhomogcneous and can be solved by standard matrix inversion or by ad-hoc methods. We empirically found 
that the numerical accuracy of the result can depend on which line is replaced by the condition on the trace. 

The Matlab® computing language has powerful built-in linear algebra and matrix manipulation tools, that allow 
translating all the above steps into a very compact code. In particular, all matrices described above are highly 
sparse. This allows using the sparse-matrix representation available in Matlab®, which results in a very moderate 
run-time memory requirement. The Kronecker product is also a built-in function in Matlab®, and takes advantage 
of the sparse matrix representation. To solve the linear system, the built-in linear solver mldivide is the fastest 
choice for small problems. It has the disadvantage however that the matrix is internally turned into its full-memory 
representation, thus rapidly hitting a memory limitation as the size of the problem grows. For larger problems, using 
one of the iterative methods available for sparse matrices is probably the best solution. We empirically found that the 
least-squares method, implemented by the Matlab® function lsqr, gives the best results for the problem treated in 
this work. In this way, the numerical solution of the stationary von Neumann problem for a given set of parameters 
required from less than one second - for the smallest problems - up to a few hours at most for the largest problems 
considered, and all computations could be carried out on a modern desktop PC with 32 GB of RAM. 



Appendix B: Weak pump limit 



In the limit of weak driving field, e — > 0, we can assume that only states with the lowest occupation numbers are 
significant for the occurrence of unconvetional photon blockade (UPB). In particular, we need to include at least 
states with one phonon and two photons, in order to account for the effective Kerr nonlinearity induced by the 
optomechanical coupling. We can then proceed in analogy with the analysis carried out by Bamba et alX in the case 
of the simple Kerr nonlinearity. We thus express the steady state of the system as 

2 1 

|V>) = X] C n in2 n m \nin 2 n m ) . (B.l) 

^1+^2—0 rim— 

To determine the twelve coefficients C nin2 n m , in the steady state, in the limit of vanishing temperature and neglecting 
pure dephasing, we can write the Schrodingcr equation for |^>) 

if=^), (B.2) 

and set d\ip)/dt = 0. Here H is an effective non-hermitian Hamiltonian that accounts for the dissipation terms, defined 
as the original Hamiltonian 

H = Aidjai + A 2 a3,a 2 + w m 5JS 2 

- J(a|a 2 + a\ai) + ga\a,2{b\ + b 2 ) + e(a[ + ai) , (B.3) 

with the replacement Aj — > Aj — in/2 and uj m — > uj m — 17/2. We further assume Ai = A 2 = A for simplicity (see 
discussion in the main article). Then, by replacing (|B.lj) into (|B.2[) . assuming steady state, and requiring that the the 
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coefficient of each number state \niri2n m ) vanish independently, we finally obtain the following set of linear equations 



= 


ACioo + eCooo + v2eC200 + JCqw , 


(B.4) 


= 


ACoio + eCuo + JCiqq + gCon , 


(B.5) 


= 


(A + w m )C n + eCin + JCioi + gCoio , 


(B.6) 


= 


(A + Lo m )C wl + eCooi + \/2eC 20 i + JCon , 


(B.7) 


= 


2AC 200 + V2eC W0 + V2JC 1W , 


(B.8) 


= 


2ACno + eCoio + V2 J(C 20 o + C 02 o) + gC ul , 


(B.9) 


= 


2AC 020 + V2JC 1W + 2gC 021 , 


(B.10) 


= 


(2A + u m )C 02 i + V2JC in + 2.gC 020 , 


(B.11) 


= 


(2A + u m )C n i + V2J(C 20 i + C ai) + gC 110 + eC 0U , 


(B-12) 


= 


(2A + co m )C 2 oi + V2eC 101 + V2JC ni , 


(B.13) 



Two additional equations, namely eCioo = and eCioi + w m Cooi = 0, are irrelevant to the problem. The first one 
vanishes identically in the limit e — > 0, while the second one determines the coefficient Cooi which however vanishes 
as 0(e 2 ). Furthermore, we can assume Cooo = 1, thus leaving the state unnormalized, which doesn't affect the value 
of g^ (0). Finally notice that the leading order in the e-dependence of the various coefficients C„ iri2 „ m is 0(e™ 1+ " 2 ). 
Then, some terms can be neglected in Eqs. (|B.4[) - (|B.13[) . as they are of higher order in e. In particular, the terms 
\/2eC 2 oo in (|B.4|) . eCuo in (|B.5[) . eCm in (|B.6[) . as well as the terms eCooi + V2eC 2 oi in (|B.7|) can be neglected. 
In this way, the matrix corresponding to the linear system is no longer hermitic as the original one, but the set of 
equations is now closed, and the coefficients C ni „ 2 „ m are given by the leading order in e only. A full analytical solution 
of the system obtained in this way is possible but rather cumbersome, and doesn't provide as much insight into the 
optimal values of parameters required for UPB, as the analogous system of equations studied in the Kerr casei Such 
an insight can instead be gathered from approximate considerations, as follows. 

The zero-delay two-photon correlation g^(0) in the first mode can be expressed as 

7 (2) rm _ (a{a\aiai) 



<7^(0) 



C 2 oo| 2 + 2|C 20 i| 2 



(|C 100 | 2 + IC101I 2 ) 2 



As all C ni n 2 n m are homogeneous in e Ill +" 2 ; this expression does not depend on e, as expected in the limit e — > 0. 
Expression (|B 141) already clarifies the occurrence of a lower bound in g( 2 \0), as discussed in the main article. It is 
determined by the presence of the one-phonon term in the numerator - which is instead absent in the corresponding 
equation for the simple Kerr model.— There is no optimal choice of the parameters that makes the two terms in the 
numerator of (|B. 14[) vanish simultaneously. In particular, the coefficient C201 never vanishes. We can prove this 
statement to lowest order in g/uj m - Let us divide Eqs. (|B.4|) - (|B.13|) by oj m and solve them iteratively in orders of 
g/uj m . Then, coefficients C niri2nm with n m = only depend on even powers of g/oj m , while those with n m = 1 depend 
on odd powers of the same parameter. More precisely, 

r A 2 g 2 J 2 (4A 2 A 1 + J 2 A 2 ) 4 

C20 ° = V2(A 2 -J 2 ) 2 + 2A(A 2 -J 2 )3(A 2 -J 2 ) + 0(5 } ' (B - 15) 
Czoi - - (A2 _ g ^_ j2) +^ 3 ) ; (MB) 

where all quantities are now normalized by w m , we have assumed 7 = for simplicity, and we have defined Ai = A + 1 
and A 2 = 2 A + 1. It is now clear that the first-order term in (|B.16[) never vanishes unless either n or J vanish. This 
is enough to prove analytically the existence of a lower bound on g^(0) in the optomechanical UPB mechanism 
occurring for g/u> m <C 1. Casting Eq. (|B.15j) in a form with a common denominator instead, we can derive the 
conditions for vanishing C 2 oo by equating to zero the numerator of such expression. This gives 

2A 3 (A 2 - J 2 )(A 2 - J 2 ) + g 2 J 2 (4A 2 Ai + J 2 A 2 ) = . (B.17) 

We can set back A — > A — in/2, Ai — > A + uj m . A 2 — > 2A + ui m , and require the real and imaginary parts of this 
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FIG. I: g (0) as a function of J computed for uj m = 24k, g = 1.16k, 7 = 0.01k, K b T = 0.1k. Line: analytical limit (|B.14J) . 
Circles: exact solution, e = 0.1k. Squares: exact solution, e = «. 



expression to vanish independently. Then, in the limit A, k, j « J « ui m , we obtain after straightforward algebra 



A op t 



^ (R18) 



= ±^^, (B.19) 

K0J m 3^3 J 

that coincide with the optimal conditions found in Ref. Q], if we assume an effective Kerr nonlinearity U c s = g 2 /u> m , 
as occurs in an optomechanical system. 

We conclude by solving the system (|B.4[) - (IB.13I) numerically (keeping only the leading-order terms in e in each 
equation as discussed above) and computing g^(0) from Eq. (|B.14[) , which holds in the limit e — > 0. Fig. [Jshows the 
quantity <7^ 2 ^(0) computed from Eq. (|B. 14|) (line), from the fully converged solution of the stationary von Neumann 
problem for e = 0.1k (circles), and e = n (squares). Clearly, the small-e result is very well reproduced by the 
approximate model presented here, which is instead unable to account for the exact result in the case of stronger 
driving field e = K. 
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